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Abstract. We assess transport properties of heavy quarks in the Quark-Gluon Plasma (QGP) that 
show a strong non-perturbative behavior A T-matrix approach based on a potential taken from 
lattice QCD hints at the presence of heavy-quark (HQ) resonant scattering with an increasing 
strength as the temperature, T, reaches the critical temperature, Tc ~ 170 MeV for deconfinement 
from above. The implementation of HQ resonance scattering along with a hadronization via quark 
coalescence under the conditions of the plasma created in heavy-ion collisions has been shown 
to correctly describe both the nuclear modification factor, Raa, and the elliptic flow, V2, of single 
electrons at RHIC and have correctly predicted the Raa of D mesons at LHC energy. 

Keywords: Quark-Gluon Plasma, Heavy-Quarks, Relativistic Heavy-Ion collisions, Quarkonia 

INTRODUCTION 

One of the most interesting questions in high-energy nuclear physics is about the prop- 
erties of the hot and dense medium created in ultra-relativistic heavy-ion collisions. 
Finite-temperature lattice-QCD (IQCD) calculations of strongly-interacting matter pre- 
dict a phase transition from hadronic matter to a quark-gluon plasma (QGP) at a critical 
temperature, Tc ~ 170 MeV [1]. The heavy charm and bottom quarks are particularly 
valuable probes for the properties of this medium. In the context of QGP physics they 
are considered heavy because their mass, mg, is large not only with respect to Aqcd but 
also to the temperature, T, of the plasma. This remains true going from SPS to LHC 
energies spanning a T range of 200-600 MeV. This property makes the study of heavy 
quark special, because mg ^ Aqcd allows to determine the initial heavy-quark (HQ) 
spectra by means of pQCD and makes available an out-of-equilibrium probe; the pro- 
duction time <^ Tqgp is much smaller than the QGP lifetime. Therefore heavy quarks 
pass through the entire evolution of the fireball; the HQ equilibration time Teq is of the 

order of the QGP lifetime but smaller than the light-quark one ~ Tqgp ^ Teq which 
means that in principle they carry more information; mg ^ T has two important impli- 
cations: on one hand it imposes effective flavor conservation, that in particular holds not 
only during the QGP phase but also in the hadronization process; on the other hand it 
implies also that the momentum exchange by collisions <^ (parametrically dom- 
inated by elastic scatterings), and the dynamics can be treated as a Brownian motion by 
means of a Fokker-Planck equation which constitutes a significant simplification of the 



study of transport properties. Finally, the three-momentum transfer dominates over en- 

ergy transfer |^| ^ ^7o ~ ^ which allows to use the concept of a potential and therefore 

to link the HQ physics to the studies to of the HQ free energy in IQCD [33]. The latter 
property allows to employ a finite-temperature T-matrix approach to study the problem 
of HQ interactions in the medium providing a consistent framework to evaluate both 
bound-state and scattering solutions based on a two-body static potential. 

The first results at RHIC further enhanced the potential interest in HQ dynamics 
by showing an unexepectedly strong interaction of heavy quarks with the medium 
observed through a small nuclear modification factor, i?^(/7j), and a large elliptic 
flow, V2{pt), of the single electrons, e^, from semileptonic HQ decays. Before these 
experimental results such a behavior was considered as an unrealistic upper limit, useful 
only as a reference [2]. Instead, the prediction of a large Raa and a small V2 obtained 
from gluon Bremsstrahlung from heavy quarks has been in striking contrast with the 
observations [4, 5]. Furthemore, it has not been possible to observe the expected mass 
hierarchy in the suppression and its color dependence that would lead to Raa{B) > 
Raa{D) > RAA{h) [4]. This was in part due to the rather indirect experimental access 
to the HQ dynamics through the measurements on the single e^, which does not allow 
to disentagle the individual contributions from B and D mesons. A first breakthorugh 
in this direction has been possible thanks to the recent preliminary results on D mesons 
presented at QM201 1 by the ALICE Collaboration [6] which confirm a large suppression 
of heavy mesons. In these Proceedings we discuss HQ scattering in the QGP focusing 
on the possibility of resonance scattering with light (anti-)quarks reminiscent of quark 
confinement. We discuss the comparison of the model with the data at RHIC and the 
2007 prediction for LHC. 



RESONANT IN-MEDIUM HEAVY-QUARK SCATTERING 

Early approaches to the study of HQ scattering in the medium was carried our using 
perturbative QCD (pQCD), first using elastic scattering [3] and later based on gluon- 
bremsstrahlung energy loss and/or including elastic HQ scattering [4, 5]. In such ap- 
proaches only a moderate decrease of Raa and a small V2 of the single electrons from 
D and B decays have been expected in clear disagreement with the experimental obser- 
vations. Hence non-perturbative approaches are expected to be necessary to explain the 
strong HQ interaction with the medium. An early suggestion postulated a mechanism 
via the formation of D- and 5-meson resonance excitations in the deconfined phase of 
QCD matter [7, 8]. This idea has first been realized through a non-relativistic effective 
field theory, modeling colorless (pseudo-) scalar and (axial-) vector D- and 5-mesons 
exploiting both chiral symmetry and HQ (spin) symmetry, 

-Gy {qt^l^-^C-q^t^,^^-^C + h.C^ 



The fields, (O*) represent (anti)-D-mesons, transforming as isospinors under isospin 
rotations and witli the usual =Sf free terms for quarks and D-mesons; see Ref. [7] for 
more details. 

The interaction terms in Eq. (1) were evaluated to leading order in l/rric according 
to HQ effective theory (HQET). The main parameter is given by the coupling Gsy 
varied to allow for widths of the D-meson spectral functions of 300-500 MeV, to 
approximately cover the range suggested by effective quark models. It is important 
to note that we assume the D-meson resonances, mo = 2.0 GeV to be located above 
the c — q mass threshold, + = 1.5 GeV, which renders them accessible in c — ^ 
scattering processes. The situation is quite different for (bound) meson states (i.e., below 
the anti-/quark threshold), where the resonant part of the scattering amplitude cannot be 
probed through c + q ^ c + q interactions (even for resonance masses close to threshold, 
thermal energies of anti-/quarks imply that the average collision energy is significantly 
above the resonance peak). 

We find that the presence of these resonances at moderate QGP temperatures sub- 
stantially accelerates the kinetic equilibration of c-quarks as compared to using per- 
turbative interactions. We have concentrated on the charm-quark case, but completely 
analogous expressions apply to the bottom sector. These approaches have been used in 
2006 to make predictions for the case of Au-i-Au at RHIC energies, and the details are 
presented in Ref. [8]. The key ingredients are that the transport cross section is appre- 
ciably larger than the pQCD ones and a hadronization mechanism that includes quark 
coalescence [12, 14, 15], leading to an enhancement of both the quark Raa and V2 toward 
a better agreement with the data [8, 16]. Other approaches have consistently pointed out 
the necessity to go beyond a simple pQCD scheme [17, 18] to account for the observed 
Raa even if these studies included coalescence especially in the < 3-4 GeV region, 
which translates to pj < 7-8 GeV at the meson level. In particular our approach appeared 
to be the only one able to describe both the small Raa and large V2{pt) simultaneuosly. 
Furthermore a prediction for Pb-l-Pb collisions at LHC energy has been presented in the 
last call for prediction proceedings [32] and will be discussed in the following. 



T-matrix Approach to Heavy-Quark Scattering 

The success of the first applications of the idea of resonant HQ scattering in the 
medium has lead to a further and more realistic assessment of the existence of resonance 
scattering. The idea of the existence of resonance scattering is indeed supported by 
IQCD on the quark correlators for both heavy and light quarks [19] as well as by Non- 
Relativsitic QCD solved on lattice for heavy quarks [20]. They both show the existence 
of a peak in the spectral function even at temperatures substantially higher than Tc 
suggesting the presence of a physical mechanism beyond a simple free scattering. 

As mentioned in the introduction, the large HQ mass allows the use of an interaction 
potential between quarks. This has the advantage that T-matrix scattering theory be- 
comes applicable, which does not rely on a perturbative expansion but is able to account 
also for moderate or even strong coupling, where resummations of large diagrams are 
necessary and realized via the standard ladder sum. An extra benefit is that one can in 



principle extract the quark potential from finite-temperature lattice QCD (IQCD), or at 
least be constrained by IQCD "data" which gives a parameter-free input. 

A static heavy-quark light-quark potential, V{r), has been used in the vacuum to 
successfully describe D- meson spectra and decays [22, 23]. We assume that the effective 
in-medium potential can be extracted from finite-temperature IQCD calculations of the 
color-singlet free energy Fi(r, T) [24, 25] for a static QQ pair as the internal potential 
energy by the usual thermodynamic relation [21, 26, 27], 

c/i(r,r)=Fi(r,r)-r^^^. (2) 

For the application as a scattering kernel in a T-matrix equation, the potential has to 
vanish for r — )■ oo. Thus we choose the accordingly subtracted internal potential energy, 

yi(r,r) = C/i(r,r)-C/i(r^oo,r). (3) 

In IQCD simulations one finds that Ui{r oo, T) is a decreasing function with tem- 
perature which could be associated as a contribution to the in-medium HQ mass, 
mg(r) = niQ + Ui{r ^ oo, T) /2, where niQ denotes the bare mass. However, close to 
Tc the asymptotic value, Ui{r ^ oo^T), develops a pronounced peak structure which is 
currently not fully understood (possibly related to multiple coupled channel effects). 
For simplicity, in the current calculation, we assume constant effective HQ masses, 
nic = 1.5 GeV and nih = 4.5 GeV. 

We also consider the complete set of color channels for the Qq (singlet and octet) and 
Qq (anti-triplet and sextet) systems, using Casimir scaling as in leading-order pQCD, 
Vs = -Vi /8, ^3 = Vi /2. , Ve = -Vi /4, which is also justified by IQCD calculations 
of the finite-r HQ free energy [28]. 

The uncertainties due to the extraction of the potential have a moderate final impact 
on the agreement with the experimental data, see also Fig. 3 (left), while the difference in 
the transport coefficients going from the internal energy, U, to the free energy, F, is quite 
substantial. However, the successful application to compute quarkonium correlators and 
HQ susceptibilities lends a-posteriori support (albeit not validation) for the choice of U, 
which is not as convincing for F (cf. Ref. [29]). 

The starting point of HQ Brueckner theory is a system of coupled Bethe-Salpeter (BS) 
and Sch winger-Dyson (SD) equations characterizing the HQ interaction and propagator 
in the QGP, 

M = K + JkGM, = 1.^ + JtS'^, = S^ + 8^11^8^, (4) 

where M denotes the scattering amplitude between a heavy (Q) and a light (q) quark 
or antiquark, K the two-body interaction kernel, G the two-particle (qQ) propagator, 

5^'^ (Sq '^) the (free) single-particle propagator, and the HQ selfenergy receiving 
contributions from thermal gluons (Eg) and light quarks (where the latter are com- 
puted self-consistently from the heavy-light scattering amplitude). Since we focus on 
a QGP at zero chemical potential (/i^=0), all quantities are quark-antiquark symmet- 
ric. The predominantly space-like momentum transfer in on-shell scattering of heavy 
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FIGURE 1. Diagrammatic representation of the Brueckner many-body scheme for the coupled system 
of the T-matrix based on IQCD static internal potential energy as the interaction kernel and the HQ 
self-energy. 



quarks, = q^ — q^ c:^ —Q^, justifies a static (potential) approximation to its interac- 
tion, K . This allows to reduce the four-dimensional (4D) BS equation into a 3D 
Lippmann-Schwinger equation for the T-matrix, which greatly simplifies its solution. 
Using azimuthal symmetry and a partial-wave expansion leads to a one-dimensional 
integral equation for the amplitudes, T/, of given angular momentum, /, 

TnE;p',p)=vnp'.p) + lJ dkk^ Vt{p',k) GQq{E-k) Ti%E-k,p) , (5) 

where we also indicated the four possible color channels, a, for Qq (singlet and octet, 
a=l and 8) and Qq (antitriplet and sextet, a = 3 and 6) states. 
The in-medium HQ quasiparticle-dispersion relation is given by 

(of = ^k2 + ml{T)+ReLQ{cof,k), (6) 

where it has been assumed that Eq. (4) (middle) can be decomposed into two distinct 
contributions: a "gluon-induced" one. Eg , generating a temperature-dependent mass, 
mg(r), and a selfenergy, E^, due to scattering of heavy quarks off thermal light quarks. 
The former is associated with the long-distance limit of the potential, while the latter 
is explicitly evaluated from the above heavy-light T-matrix. To close the system of 
equations (4) in the quark sector, the analogous system for the light sector is required, 
which has been solved selfconsistently for E^ and Tqq^qq in Ref. [21]. Guided by the 
results obtained there we employ a light-quark propagator (figuring into the second 
term in Eq. (4) with a constant thermal light-quark mass, niq — 0.25 GeV, and width, 
r^ = 200MeV = -2E'?. 

Clearly, the potential approximation is less reliable in the light-quark sector, but it 
turns out that their in-medium self energies (real and imaginary parts), which are needed 
for the HQ selfenergy, Eq. (4), have a small effect on both the scattering amplitude, TQq, 
and the HQ transport coeffcients. 

The above system of equations (for T and V) is pictorially represented in Fig. 1 with 
the upper, middle and lower panel corresponding to Eqs. (4), respectively. 
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FIGURE 2. Left: Imaginary part of the S-v/ave in-medium T matrix for cq and cq scattering in the 
color-singlet and -antitriplet channels based on the parameterization of the IQCD potential energy by 
[Wo]. Right: The drag coefficient, 7, as a function of HQ momentum, calculated via (4) with scattering- 
matrix elements from the non-perturbative T-matrix calculation (using the parameterization of the IQCD 
internal potential energies by [Wo]) compared to a LO perturbative calculation based on matrix elements. 



We restrict ourselves to 5 (/ = 0) and P (I = 1) waves. As can be seen from Fig. 2 
(left), in the dominating attractive color-singlet Qq and color-antitriplet Qq channels, 
close to the critical temperature, Tc, resonance states close to threshold, E'thr = i^q + mq 
are formed, similar as conjectured in the effective resonance model described in the 
previous section [7, 8]. However, in this full in-medium scheme the resonances melt at 
lower temperatures T > IJ Tc and T >lATc, respectively. 



HEAVY-QUARK OBSERVABLES IN HEAVY-ION COLLISIONS 

A direct comparison between the microscopic description of the HQ dynamics in the 
QGP and the experimental observables necessitates a dynamical implementation of 
the HQ scattering in the medium plus a model for the hadronization and finally the 
implementation of the semileptonic decay. 

The HQ motion in the hot and dense QGP, consisting of light quarks and gluons, can 
be described by a Langevin simulation of the Fokker-Planck equation. 

The drag or friction coefficient, 7, and diffusion coefficients, Bij, are calculated from 
the invariant scattering-matrix elements [3]. Taking into account elastic scattering of the 
heavy quark with a light quark or antiquark given in terms of the above calculated T- 
matrix one can calculate the drag with a standard procedures the HQ drag and diffusion 
coefficients, shown in Fig. 2, see Ref. [34] for more details. 

The nonperturbative HQ light-quark scattering-matrix elements are supplemented by 
the corresponding perturbative elastic HQ gluon- scattering ones. The ^-channel singu- 
larity is regulated by a gluon-Debye screening mass of nig = gT with a strong coupling 
constant, g = y/Anas, using = 0.4. 



As shown in Fig. 2, close to Tc the equilibration times, Teq = I/7 ~ 7 fm/c, for charm 
quarks are a factor of ~ 4 smaller than for a corresponding pQCD calculation, reminis- 
cent to the results based on the model assuming Z)-meson like resonance states above 
Tc [7, 8]. In contrast to this and other calculations of the HQ transport coefficients, here 
the drag coefficients decrease with increasing temperature because of the "melting" of 
the dynamically generated resonances at increasing temperatures due to the diminishing 
interaction strength from the IQCD potentials. 

To solve the Fokker-Planck equation (7) under conditions of the sQGP medium 
produced in heavy-ion collisions, we use an isentropically expanding thermal fireball 
model, assuming an ideal-gas equation of state of Nf — 2.5 effective massless light-quark 
flavors and gluons. The initial spatial distribution of HQ production is determined with 
a Glauber model. The c-quark spectra are taken from a modified PYTHIA calculation 
to fit D and D* spectra in d-Au collisions [35], assuming 5-function fragmentation. 
The b-quavk pr spectrum is taken from PYTHIA assuming a cross-section ratio of 
a^^/ Occ — 5 • 10~^ and a crossing of the c- and ^-decay electron spectra at pt ^5 GeV, 
consistent with FONLL pQCD calculations [31]. 

The last step toward a comparison of the above described model for HQ diffusion 
in the QGP with the data from RHIC is the hadronization of the HQ spectra to 
D- and B-mesons and their subsequent semileptonic decay. Here we use the quark- 
coalescence model described in [12, 14]. In recent years, the coalescence of quarks in 
the hot and dense medium created in heavy-ion collisions has been shown to provide 
a successful hadronization mechanism to explain phenomena such as the scaling of 
hadronic elliptic-flow parameters, V2, with the number of constituent quarks, vj^hiPt) — 
nhV2,qiPt /fih), where n;, = 2(3) for mesons (hadrons) denotes the number of constituent 
quarks contained in the hadron, h, and the large p/n ratio in Au-Au compared to p-p 
collisions [12, 13, 14]. Quark coalescence is most efficient in the low-pr regime where 
most c and b quarks combine into D and B mesons, respectively. To conserve the total 
HQ number, we assume that the remaining heavy quarks hadronize via (5-function) 
fragmentation. 

As shown in Fig. 3 the Langevin simulation of the HQ diffusion, followed by the 
combined quark-coalescence fragmentation description of hadronization to D and B 
mesons and their subsequent semileptonic decay, successfully accounts for both the Raa 
and V2 of in 200AGeV Au-Au collisions [11, 10] at RHIC. The uncertainty due to 
two different parameterizations of the IQCD potentials by [Wo] [27] and [SZ] [26] is not 
so large. 

Comparing the solid and dashed lines one sees that the effect from the "momentum 

kick" of the light quarks in coalescence, an enhancement of both Raa and V2, is important 
for the agreement of both observables with the data. As can be seen from the lower panel 
in Fig. 3 (right), the effects of the 5-meson decay contribution to the spectra become 
visible for pr > 2.5-3 GeV. A closer inspection of the time evolution of the pt spectra 
shows that the suppression of high-pj- heavy quarks occurs mostly in the beginning 
of the time evolution, while the V2 is built up later at temperatures close to which 
is to be expected since the V2 of the bulk medium is fully developed at later stages 
only [37, 38]. This effect is more pronounced due to resonance formation, because the 
transport coefficients become larger close to Tc, or at least decrease much slower than in 
the pQCD case [29]. 
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FIGURE 3. Left: Results for the nuclear modification factor (upper panel) and elliptic flow (lower 
panel) of single electrons with/without (solid/dashed lines) quark coalescence in Au-Au collisions com- 
pared to RHIC data [10, 11]. Right: As in the left panel but only for the case of [Wo] potential with 
coalescence included it is shown the contribution of single e± from D and B decay. 



Prediction at LHC 

In 2007 we made predictions for the the D and B spectra at the LHC employing the 
effective resonant model described above. For D mesons we know that the model is quite 
reliable in the sense that it generates predictions similar to the T-matrix approach. Since 
the initial temperatures at the LHC are expected to exceed the resonance dissociation 
temperatures, the prediction implemented a switching off of the resonances at r^iss = 
ITc = 360 MeV by a factor (1 +exp[(r - Tai,,)/A])-^ with (A = 50 MeV) in the 
transport coefficients. 

The temperature evolution in the fireball assumed a total entropy fixed by the num- 
ber of charged hadrons which we have extrapolated to dNch/^y — 1400 for central 
■s/snn = 5.5 TeV Pb-Pb collisions leading to an initial temperature, Tq ~ 520 IMeV. This 
multiplicity truns out to be very close to the one measured at y^s^yyv = 2.75 TeV, hence 
now we know that these have to be considered more properly as predictions for this 
energy. 

Initial HQ pr spectra are computed using PYTHIA with parameters as used by the 
ALICE Collaboration [36]. Hadronization is treated as previously discussed for the 
RHIC case. The shadowing has not been included but as shown by the orange solid 
line this should not affect the prediction at pj- > 3.5 GeV. 

In Fig. 4 the predictions for the D-meson Raa are shown by the shaded red area 
corresponding to the uncertainties in the resonance model. The data from the ALICE 
Collaboration [6] are shown by rectangles, red for D+ and green for D^, and are for 
the — 20% centrality and therefore corresponding to an impact parameter, b = 4.7 fm, 
hence more central with respect to the 2007 calculations. For this reason in the figure 
the red dashed line is drawn to indicate the extrapolation from the b = 1 fm calculation 
to the more central with b = 4.7 fm which results in a 15% correction. It is clear that 
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FIGURE 4. The nuclear modification for D meson at LHC for Ph + Pb at ^^snn = 2.75 TeV for more 
detail see the text. 

the prediction from simple pQCD elastic scattering, shown by the blue dash-dotted line, 
cannot account at all for the observed suppression of the D spectra. 

Data with a better statistics and especially a measurement of the elliptic flow will 
greatly improve the discrimination power of the model. From the theorethical side it is 
desireable to have predictions using the T-matrix approach based on the IQCD-based 
potential, especially for the B mesons. These are differently affected by the medium and 
should show a resonant scattering persistent up to higher temperatures compared to the 
D mesons resulting in a smaller Raa compared to the simplified resonance model shown 
by the violet shaded area shown in Fig. 4. 

CONCLUSIONS 

The in-medium interaction of heavy quarks have been intriguing from the first prelim- 
inary results at RHIC showing a much larger suppression than expected, together with 
a large elliptic flow, as inferred from semileptonic HQ decay electrons. An in-medium 
T-matrix approach (utilizing potentials estimated from lattice QCD) has been applied to 
evaluate HQ interactions in the QGP showing the existence of in medium prehadronic 
and diquark resonance states, which increase in strength when approaching Tc. When 
implemented into Langevin simulations at RHIC, reasonable agreement with both the 
suppression and elliptic flow of spectra from HQ decays emerges. 

Several theoretical problems remain open like a proper definition of the potential, 
corrections to the T-matrix approach including radiative ones, in-medium mass and 
width effects as well as a self-consistent treatment of the hadronization via coalescence. 
On the other hand this field is entering a new stage thanks to the possibility to disentangle 
the B and D mesons which in itself can provide the key to discriminate different models 
[16, 17, 18, 39, 40] . 
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